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Abstract 

The collisions of cosmic strings loops and the dynamics of junctions formations in expanding 
backgrounds are studied. The key parameter controlling the dynamics of junctions formation, 
the cosmic strings zipping and unzipping is the relative size of the loops compared to the Hubble 
expansion rate at the time of collision. We study analytically and numerically these processes for 
large super-horizon size loops, for small sub-horizon size loops as well as for loops with the radii 
comparable to the Hubble expansion rate at the time of collision. 
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I. INTRODUCTION 



In models of brane inflation cosmic strings are copiously produced [l, 2|, for reviews see 
e.g. ^-7\. These cosmic superstrings are in the forms of fundamental strings (F-strings), 
Dl-branes (D-strings) or the bound states of p F-strings and q D-strings, the (p, q) strings. 
When two {p, q) cosmic superstrings collide junctions are formed due to charge conservation. 
This is in contrast to the collision of conventional gauge strings where upon collision they 
exchange partners and intercommute with the probability close to unity. Therefore one may 
consider the junction formation as a novel feature of a network of cosmic superstrings which 
may prove crucial in cosmic superstrings detection in cosmological observations. Networks 
of cosmic strings with junctions have interesting physical properties, such as the formation 
of multiple images and non-trivial gravitational wave emission 
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Different 



theoretical aspects of (p, q) string construction were studied in |12|-|l7l| while the cosmological 



18|. 



evolution of a string network with junctions has been investigated in 

In a recent paper [l9| the collision of two loops of cosmic strings in a flat background 
was studied. It was found that with appropriate initial conditions determined by the angle 
of collision, the colliding loops velocities and the loops relative tensions, junctions can form. 
However, after the junction is formed it can not grow indeflnitely and after some time the 
junction start to unzip and the colliding loops disentangle and pass by from each other. The 
junctions' zipping and unzipping are interesting and yet non-trivial dynamical properties. 
These phenomena becomes more signiflcant in the light of cosmic strings simulation by 



Urrestilla and Vilenkin 



20|. In their model, the cosmic strings are two types of U(l) gauge 



strings with interactions between them. Due to the interaction, the strings cannot exchange 
partners and a bound state will form if the strings are not moving too fast. It was shown 
that the length and the distribution of the string network are dominated by the original 
strings and there is a negligible contribution to the string network length and population 
from the bound states strings. This can be understood based on the following two reasons. 
Firstly, the junctions may not form if the colliding strings are moving very fast so they can 



simply pass through each other 
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Secondly and more interestingly, if the junctions 



are formed, they start to unzip during the evolution. 



Our aim here is to generalize the results of 



19j to the case of cosmic strings loops collision 



in cosmological backgrounds, i.e. the radiation and the matter dominated era. As we shall 
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see in sections UTIl and HVl the size of the loops compared to the Hubble radius at the time 
of collision plays a significant role in junctions evolutions and cosmic strings zipping and 
unzipping. 

The paper is organized as follows. In section [TTl we present our set up and provide the 
formalism of junction formation for arbitrary cosmic strings loops colliding in cosmological 



backgrounds. This is a generalization of [22| and [28| where they presented the formalism of 
cosmic strings collision in the flat background. In section [XTTl we concentrate to the example 
of two identical loops in cosmological backgrounds. After setting the background equations 
for the loops profiles, we present the equations governing the dynamics of the junctions. 
For the case of very large loops (super-horizon size loops) and very small loops (sub-horizon 
size loops) we are able to present some analytical results. In section IIVI we present our full 
numerical results for different loop configurations. The conclusion is given in section |V] 

II. LOOPS COLLISION IN AN EXPANDING UNIVERSE 

Here we present the formalism of junction formation for arbitrary cosmic strings loops 
colliding in a cosmological background. In section IIIII we employ the results obtained in 
this section to the particular example of two identical loops at collision in cosmological 
backgrounds. The formalism of cosmic strings collision in a flat background was studied in 



22| and |28|. 



Our cosmological background is the standard FRW metric 

ds^ = a^{T){dr^ -dx"), (1) 

where r is the conformal time related to the cosmic time t via dt = adr, a(r) is the scale factor 
and we assume that the background space-time has no spatial curvature. Our cosmological 
background is either radiation dominated (RD) or matter dominated (MD). 

Suppose represents the profile of the i-th cosmic string in the target space-time. As 
usual, we can go to the temporal gauge where the time on the string world-sheet is the same 
as the conformal time, X° = r, and Xf = (r, Xj). Denoting the other coordinate of the 
world-sheet parameterization by a, the gauge condition 

x,-x: = 0, (2) 

holds where the ■ and the prime indicate the derivatives with respect to r and a respectively. 



A schematic view of two loops of cosmic strings in collision is shown in Fig. [T] . After 
collision, there are four junctions and eight kinks. The formation of the kinks is a manifes- 
tation of the fact that the speed of light is finite and parts of the old strings which did not 
"feel" the formation of junctions evolve as before. In the following we denote the incoming 
strings by x, where i = 1,2 whereas the newly formed strings are denoted by where 
a = 1,2,3. The junctions and the kinks on each string are described by cTq = Sai^r) and 
(Ji = Wj(r) respectively. 



As mentioned in 
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one complexity of dealing with loops in collision is the orienta- 



tion of the (jj coordinate at junctions. We follow the prescription of |28| and use the sign 
parameterization for 6^ according to which 6^ can take values ±1. If the value of cTa of a 
particular string increases(decrease) towards the junction J, we assign 6;^ = +1(5^ = —1). 
With this prescription, the two ends of a piece of string ending in two neighboring junctions 
have opposite 6 parameters. The arrows in Fig. [1] indicate this prescription. Since it is 
important for the later analysis, we now give the values of 6^ at each junction: 





6i = +1 




5i = -l 




61 = 


-1 






+1 




A : 


§2 = -1 


B : 


S2 = +1 


C : 


62 = 


+1 


D : 


62 = 


-1 


(3) 




53 = -1 




S3 = +l 




S3 = 


+1 




53 = 


-1 





Equipped with the ^-prescription, the action for the system of two loops in collision is 
S = - ^ X^/i J dTda,a\T)^{l-yf)y'M6^{si{r) - a,)) 9(6^ {a. - a;/(r))) (4) 

J i=l 

i=i 

- X^/i. / drda,a\r)^J{l - xf) ^,'e{5? {uj't {r) - a,)) e{5^{uj^{r) - a,)) 
i=i ^ 

- /i3 / dTda,a'{T)^{l - yi) y'^^iisiir) - a,)) e{6f{s^{T) - a,)) 



- /xay drdasa\T)^il - yf) y^^^(53^(.f (r) - ^3)) ei6^{s^{T) - a^)) 
+ EE / dTa\T)i^.[y^isi{T),r)-rir)] 

J a=l 

+ EE / ^™^(r)kf ■[x,(a;/(r),r)-y,(u;/(r),r)] 



J 1=1 



FIG. 1: A schematic view of the loops at the time of conision(left) and after colhsion (right). The 
arrows in the right figure indicate the directions in which the ai coordinate increases. We use the 
convention that on a loop ai runs counter clockwise. There are four junctions and eight kinks in 
total. 

where J represents the junctions A, B, C and D collectively. Here f/ and are Lagrange 
multipliers which enforce that at the kinks the newly formed strings and the old strings 
meet, Xj(u;/(r),r) = yj((X'/(r), r) and on the junctions the three newly formed strings join 
together ya(s;((r), r) = y"'(T) where y'^(r) represents the position of the junction J in target 
space. As described above, in our convention 5/ is +1 if a increases towards the junction 
and -1 in the opposite case. 

Varying the action with respect to f/, k;[ and y;^, respectively, results in 



Varying the action with respect to Xj results in the following standard equations [29] for 
the segments of old strings extended between two nearby kinks which are not influenced by 
the junctions formations 




(5) 
(6) 
(7) 





a 



whereas matching the Dirac delta functions gives the following boundary conditions at the 
kinks (jj = uji {t) 

kf = /i, ('i.e.A^^/ + , cr. = a;/(r). (9) 



Here is defined by [29 1 



x'2 



(10) 



Similarly, varying the action with respect to results in the following equations for the 
newly formed strings stretched between a junction and the nearby kink 

|(y..,J + 2^y..,.^^(^). (U) 

Now there are Dirac delta functions at the junctions da = s'^{t) and the kinks (for strings 
1 and 2) cr.j = which result in the following boundary conditions 

kf = /i. (^y.e^A''^/ + ^5/^ , a. = u;/(r) (12) 

i^ = f^a(yaeyjiSi + ^6A , aa = siir). (13) 

Combining Eqs and f lT5]) result in 

E-^A'fyaeyjf+f^) =0. (14) 

Eliminating k/ from Eqs. ^ and (IT^ . and using the gauge condition ([2]), it is easy to 
show that 

ex, = ey, = (15) 
bl^ie, = -1 (16) 

The solutions of the loops in cosmological backgrounds can not be expressed in terms of 
the the usual right- and left-movers. However, one can define the right- and left-momenta 
Pa as 

ya , -^"aJa 



p^i = - ± (17) 
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with p^^ = 1. Starting with the time derivative of ([5]) 

x'w/ + = y'ui + y, , (18) 

one can show that p~y. = P~Xi- This is the key formula which relates the unknown quantities 
Py^"^ for the newly formed strings to the known quantities Px^"' from the old strings. 

Starting with the time-derivative of Eq. (|5]) combined with Eq. f lT4|) one obtains an 
equation for Py^ in terms of Py^^ 

5i{l + 5ieJi) p+ = - 5ieX) P" - - - 5ie,si) p;^ . (19) 

^ h 

Imposing the condition p^^ = 1 one obtains 

9a - ga ^ P'bQb Cab + ^ fibficQbgc Cfoc - 1 = (20) 
6 be 

where fia = f^Jfi, = T^af^a, qi = I - ^i^asi and Cab = ^i^iVy, ■ Py,- 

Now we provide an equation for the energy conversation at the junction. For this purpose, 
multiplying Eq ( !20l) by Jia and summing over a results in 

P'aQa - ^ Jj^ajj'b Cab + ^ P^afibficQbqc Cafe - 1 = (21) 
a ab abc 

However, with the reshuffling of the indices, one can easily check that the second and the 
third term above cancel out and one obtains 



or 



^/^aga = l (22) 



^5i^laeaSa = ^. (23) 



This is a generalization of the case of strings in flat background studied in [22] corresponding 
to da = 0. 

Deflning Cab = ^ — Cab and using the energy conservation one can check that Eq. ( I20ll 
results in 

qa ^ fibCabqb = ^ fibP'cCbcqbqc (24) 

fe be 
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Writing explicitly, and using ( 123|1 . yields 



93(1 



^2(1 



2/^3)(Ai2g2C23 + /^igiCls) 



2/i2)(/iigiCi2 +/i393C23) 



2/x2Aiig2gici2 . 



2/U2/i3g2g3C23 



2/ii/i3gig3Ci3 



(25) 
(26) 
(27) 



Eliminating 53 and ^2 from the first two equations above and plugging in fl23|) results in our 
main formula of interest 



where Mi = fif — [^2 — f^s)"^ with a similar definition for M2 and M3. One can also obtains a 
similar equation for ^2,3 with an appropriate permutation of the indices. This set of equations 
for s'^ is our starting point to study the evolutions of junctions. 

III. JUNCTION EVOLUTIONS 

In previous section we have presented the general formalism of cosmic strings loops col- 
lision in a cosmological background. Here we specialize to the example of two identical 
loops at collision in a cosmological background where the analysis can be handled somewhat 
analytically. 

Before dealing with the loops in collision, here we summarize the background solutions 
for a loop in expanding background. Suppose the collision happens at the time r = tq. One 
can check that the cosmic time t and the conformal time r are related by t = r"+^/r|^(?2+ 1) 
where we have considered a power law expansion for the scale factor a(r) = (t/tq)". For 
a radiation an matter dominated universe, n = 1,2 respectively. In this convention, the 
scale factor at the time of collision is equal to unity. Also, calculation the Hubble expansion 
rate, H = da/adt, one can check that the Hubble expansion rate at the time of collision is 
^0 = n/TQ. 

Consider a loop extended in x — y plane moving relativistically in z direction. We choose 
the following ansatz for the loop configuration 



1 - 



J 

1 



/iMiC23 



(28) 



/ii [M1C23 + M2C13 + M3C12] ' 




(29) 
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In this picture, -R(t) = a{T)f{T) is the physical radius of the loop and Rq = /(tq) represents 
the size of the loop at the time of collision. 
The independent equations of motions are 

F" + -v^- + (1 - _ F'^)F-' = (30) 

X 

Here the loop center of mass velocity is defined by f = z{t). For the ease of the numerical 
investigations, we introduced the dimensionless time variable x = t/tq and F{x) = f{T)/TQ. 
Also the prime here and below represents derivatives with respect to the dimensionless time 

X. 

This definition leads to F = fHo/n which has a simple physical interpretation as follows. 
In our convention a{x) = x", so at the time of collision, corresponding to x = 1, a = 1. The 
physical radius of the loop at the time of collision therefore is Ro = f{x = 1). Therefore, 
the initial condition F{x = 1) = f{x = l)Ho/n is a measure of the physical radius of the 
loop compared to the Hubble radius at the time of collision. For Loops of super-horizon size 
at the time of collision F{x = 1) > 1, whereas for small sub-horizon sized loops at the time 
of collision F{x = 1) < 1. In the following, to simplify the notation we set F{x = 1) = Fi. 

In general it is not easy to find analytical solutions for the set of equations (1301) and (I3T1) . 
We solve this equation numerically. Our goal is to calculate Cab for the loop with ansatz (!29|) 
and then obtain the evolution of junctions However, one may get some useful analytical 
information in some certain limits. These include very large super-horizon size loops at the 
time of collision [Fi ^ 1) and small sub- horizon size loops {Fi <^ 1). 



The collisions of loops in Minkowski background was studied in 191]. Here we generalize 
that study to the case of strings loops collision in an expanding background. In order to 
simplify the analysis, we consider the symmetric case where the two incoming loops have 
equal tension and physical radius at the time of collision. We assume that the loops are 
extended in x — y plane and are moving along the z-direction with velocity ±v. A schematic 
view of this example is given in Fig. [T] . By symmetry the newly formed string 3 will be 
static extended either along x or y directions. Whether it is a x-link or a y-link junction 



depends on the angle of collision a 22]. For small enough angle of collision it is a y-link 
while for a large enough angle of collision it would be an x-link junction. To be be specific, 
we consider a y-link junction where the string 3 is extended along the y-direction. 



The profiles of tlie colliding loops are given by 



^T& + /(r)cos|^^ 



±z{t) 



(32) 



where for i=l (i=2) we choose upper sign (lower sign). Here 26 is the impact factor. 

From the continuity of the left moving momenta one has Py/ = Px^^ which can be served 



to find p . From f[T7j) we have 



( -VI - F'"^ - ^2 sin ^ - cos ^ \ 



_ J _ J 



Using c;^^ = 1 — b'j^bl^'yj " Py/ at each junction and the fact that 5/ 
obtains 



cos 



^ - sin ^ 



v 



T5/w(r) 



(33) 



-62 = — ^3 , one 



'11 



-12 



-13 



--23 



-22 



-33 







l + v^-{l- 2F'2 - v^) cos 2^/ + 2S(FWl - F'^ - t;2 sin 2^/ 
1 + Vl - F'2 - t;2 COS S'i - 6(F' sin 



-23 



+ C13 



1 - Vl - F'2 - COS - 6(F' sin 

2 + 2VI - F'2 - cos 5/ - 26(F' sin 



(34) 
(35) 
(36) 
(37) 
(38) 



Here and below, to simplify the notation the definition 5";^ = s'^/Rq is introduced. We 
note that 5*;^ is dimensionless which is more suitable for numerical analysis. To obtain 
Eqs. (EMSZl) we note that due to symmetry in problem, one has S'/ = —S'2 at each 



junction J. On the other hand, one also observes that Sf{x = 1) + 6*2 (x = 1) = vr and 



Si{x = 1) + 5*^(0; = 1) = Stt which was used to simplify the final results in Eqs. fl35l l3 
Plugging Cab in our master equation (128|) . the evolution of the junction is given by 



^3 {k - l)(ci3 + C23) - KC12 



Fi (k - l)(ci3 + C23) + C12 

where the dimensionless parameter k is given by the ratio of the tensions k = 2^i/ 
Also from the energy conservation Eq. ( 123|) one has 

Fi 



(39) 



S'{ 



kF 



Vl - F'2 - S'' . 



(40) 
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Equations (l39ll and (HOll jointly can be used to solve for 5*3 . Due to symmetries involved 
in the problem, we only need to find the evolution of junctions B and D and the evolutions 
of junctions A and C are mirror images of junctions B and D. We note that the conditions 
for the junction formation is that the string 3 stretching between junctions B and D to 
be created. This requires that its length to increase initially with time: S'^ > where 
measures the length of string BD. As we shall see in our numerical results, 
usually this is translated into S^' > and S^' < 0. However, in some very fine-tuned 
situations one can also find examples where S^' > and S^' > such that 5*3 > is still 
satisfied. 

I — I 

As explained in pj;!], after the junction formation, the entangled loops start to unzip. 
The onset of unzipping at junction J happens when 5*3'^ vanishes and changes its sign. As 
we shall see later, the unzipping times for junctions B and D are not equal. Since we are 
mainly interested in the evolution of the newly formed string BD, we define the onset of 
unzipping for string BD when 5*3 reaches a maximum and 5*3 = 0. After that the length of 
string BD reduces with time. Sometime after unzipping, the loops disentangle from each 
other and pass by in opposite directions. The time of loops disentanglement happens when 
the junctions B and D meet corresponding to 5*3 = 0. However, we also encounter examples 
where the loops shrink to zero before they disentangle from each other. 

As explained above, the onset of unzipping at junction J is determined when 5*3 = 0. 
Here we show that the denominator in Eq. fl39|) is always positive so the sign of 5*3'^ is 
controlled by the numerator of the above expressions. To see that note that — ^ 
pcosO + gsin^ < \/j^~+~q^ for real numbers p and q and arbitrary angle 6. Using these 
inequalities, one can easily check that C12 > 2v^ and C13 + £23 > 2(1 - Vl-v^). On the 
other hand, as described in {22!, one also requires that > ^3 for the junction formation to 
be allowed kinematically. In conclusion the denominator in 5*3"^ expression is always positive 
and the sign of S'^"^ evolution is determined by the numerator of the above expressions. 
This plays important rules in determining the junctions unzipping times in the following 
discussions. 

As explained before, our goal is to solve the background loop equations ( l30i) and ( 13T1) and 
use the resulting values of F{x) and v{x) in Cab expressions to find the junctions evolutions 
from Eq. (l39|l . This procedure can be done only numerically because both F{x) and v{x) 
can not be found analytically in general. Before presenting our full numerical analysis, we 
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consider two different limits where some analytical insights can be obtained for the junctions 
evolutions. 



A. Large super-horizon size loops 

Here we consider the limit where the colliding loops are much larger than the Hubble 
radius at the time of collision, Fi ^ 1. For the super-horizon size loops, one expect that they 
are conformally stretched as the universe expands, R{t) oc a(r), where F' and F" are small 
compared to unity. We will demonstrate this in our numerical analysis. Due to damping 
effects from the expanding background, the super-horizon loops become non-relativistic and 
V <^ 1. In this approximation, one can easily solve the background equations ( 130|) and ( 13T|) 
. Denoting F = Fi + A where A represents the small evolution of F, one obtains 

A" + -A' + ^^0, 
X Fi 

which after neglecting the sub-dominant term, results in 

2 

A ~ F' ~ - (411 

2(1 + 2n)Fi ' (1 + 2n)Fi ' ^ ^ 

As the universe expands, the loop reenter the horizon in its subsequent evolutions. This 
can be approximated when A ~ —Fi which results in the time of the loop horizon reentry 

X Ik 



^2{l + 2n)Fi . (42) 



Our full numerical analysis, as we shall see in next section, verify that this is indeed a good 
approximation. From this expression for we see that with similar initial conditions, it 
takes longer for the loops to reenter the horizon in a matter dominated era as compared to 
the radiation dominated era. 

As the super-horizon size loops stretches conformally, its center of mass velocity reduces 
rapidly. For time smaller than x^, one can find an approximate solution for v{x). Neglecting 
the terms containing f ^ and F'"^ in Eq. ( |3Tll one obtains v' + 2nv/x ~ which easily can 
be solved to give 

- ^ ' (43) 
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where vi is the value of v{x) at the time of colhsion, vi = v{x = 1). As explained above, 
we see that the loop central mass velocity reduces rapidly with time. We also see that for 
the matter dominated background with n = 2, the loss of velocity is more pronounced as 
compared to the radiation dominated background with n = 1. 

For the super-horizon size loops and keeping only terms up to F' in Cab and neglecting 
F'"^ and f ^ as explained above, one obtains 

ci2 ~ 1 - cos 2S( + 2S(F' sin 25/ , £13 + C23 ^ 2 + 2 cos S( - 25(F' sin S( 

Plugging these into S'^'^ expression results in 



, J ^ (1 + cos 5/) {k cos 5/ - 1) - 5(F' sin 5/ {k - 1 + 2k cos 5/) 

~Yi{l + cos 5/) (- cos + k)- 6{F' sin (k - 1 - 2 cos 5/) ' ^^^^ 

At the time of collision, one can neglect the terms containing F' in above expression and for 

small X one obtains 

^' R / -I \ ^ COS Si — 1 / J -\ 

X ~ 1 ~ ^—j- . (45) 

For the junction to form we need S^'{x = 1) > 0. With S'f (x = 1) = a, the condition for 
junction formation is translated into 

< a < ac , Oic = arccos(fi;~^) . (46) 



Interestingly, this is identical to the bound obtained in [23| for straight strings in collision at 
the small velocity approximation. This is expected, since the super-horizon size loops can 
be locally well approximated by the straight strings. Our full numerical analysis, presented 
in next section, indeed show that the bound on Uc given by the above equation works very 
accurately for the large loops. 

As time goes by, the terms containing F' in Eq. (jH]) becomes important. We note that 
the denominator in Eq. fH4|) is positive so the sign of 5*3 evolution is determined by the 
numerator of Eq. (jH]). Consider junction B for example. For the junction B to unzip, 
should slow down, requiring that the term containing F' in numerator of Eq. (1441) gives a 
negative contribution. With 5f = — 1 and F' < the F' correction in numerator of Eq. (jHj) 
indeed contributes negatively. This indicates that as time goes by, the rate of evolution of 
slows down until S'^^ = when the junction B starts to unzip. Similar argument applies 
to junction D too. 
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B. Small sub-horizon size loops 



Now we consider the limit where the colhding loops are much smaller than the Hubble 
radius at the time of collision, Fi <^ 1. In this limit, the damping terms in Eqs. fISU]) can 
be neglected 29(] and the loop evolution is the same as in the flat background. In this limit 
V is nearly constant and the loop has a simple periodic profile 

' X — V 



F{x) ~ Fi cos 



7F1 



(47) 



where 7 = 1/ ^/l — v"^ is the Lorentz factor. To simplify the analysis, here we chose the initial 
configuration such that F' = at the time of collision. Neglecting the effects of expansion, 
one would expect the criteria for zipping, unzipping and the loops disentanglement would 
be similar to cosmic strings loops collision in a flat background studied in |l9|. 
Starting from the energy conservation formula ( l40l) . one obtains 



St 



B{D) 



s. 



B{D) 



+ a . 



(48) 



Also calculating Cab yields 



C12 



27 



-2 2 
COS 



S( + 6., 



,x — \ 



C12 



+ = 2 + 27"^ COS ( S( + 6. 



r X 



' 7^1 



(49) 



Plugging these in S'^"^ expression yields 



S: 



I J 



s:J KCOS 
"3 



— + + a — 7 



K7 — cos 



K7 7F1 ' ^ 



(50) 



The details of the loops zipping, unzipping and disentanglement were studied in 19|]. Here 
we briefly outline the main results. For the junctions to form one requires that < a < ac 
where ac = arccos(7/fi;). The unzipping times for junctions B and D, x^ and x^, satisfy 



•Ji' •Ji' 



B 



27aFi (1 - K-2) ^ 



1 sin a 

/t7 a 



(51) 



Since sin a/a and I/K7 are always less than unity, one concludes that > x^, indicating 
that the junction B which holds the external large arcs unzip sooner than the junction D 
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X 



FIG. 2: Here the background evolution of Eqs. (f30]l and ([31]) are presented with Fi = 100, vi = 0.2 
and F'{x = 1) = 0.1. The left figure shows F{x) whereas the right figure represents v{x). The 
solid lines are for the radiation dominated backgrounds (n = 1) and the dashed lines are for the 
matter dominated backgrounds (n = 2). We are considering the loops evolution until they shrink, 
i.e. until the first root of F{x) = 0. 



which holds the internal smaller arcs. Although we have proved this only for small sub- 
horizon loops but our numerical analysis show that this conclusion holds true in general. 

The loops disentangle at the time Xf when 5*3 = 0, which is given by the parametric 
relation 



K'J COS 



cos 



Xf 



where 



Vl — — K7 /iia + sin a = , 



{Xf 



(52) 



kFi sin 



This is an implicit equation for Xf which should be solved in terms of k, 7, a and Fi. For 
this to make sense, we demand that Xf — 1 < 7Fi7r/2 before the loops shrink to zero. 



IV. NUMERICAL ANALYSIS 



In this section we present our full numerical results for different loops configurations. To 
be specific, we consider three examples of (a): large super-horizon size loops with Fi = 100, 
(b): intermediate size loops with Fi = 0.5 and (c): small sub-horizon size loops with 
Fi = 0.01. 
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FIG. 3: In these plots we have presented the evolution of junctions B and D for Fi = 100, a = tt/9 
and K = 1.2. In the left figure the upper solid red curve represents whereas the lower solid 
blue curve is that of for the radiation dominated background. The dashed curve represents 
the corresponding curves in the matter dominated background. The right graph represents the 
length of the newly formed strings fi-^, S-^ = — S^. The solid (dashed) curve is for the radiation 
(matter) dominated background. We see that in both backgrounds, after junction formation, the 
string BD reaches a maximum length and get unzipped. However, only in the radiation dominated 
background S3 becomes zero before loops shrink indicating the loops disentanglement. 

A. Large super-horizon size loops 

For super-horizon size loops with Fi ^ 1 one expects that the loops are conformally 
stretched until they re-enter the horizon. In this period, the loops lose much of its center 
of mass velocity as demonstrated by Eq. fHSl) . As mentioned before, the loss of velocity is 
more significant for the matter dominated backgrounds. Also from Eq. (H2l) we see that 
it takes longer for the loop to shrink for the matter dominated backgrounds as compared 
to the radiation dominated backgrounds. Both of these analytical conclusions were verified 
in our full numerical investigations. In Figj2] we have presented the background solutions 
of F{x) and v{x) solving Eqs. fl30l) and fl3T|) numerically. As is clear from the figure, in a 
matter dominated background, it takes longer for the loop to shrink. This in turn plays some 
roles in the junctions evolutions and loops disentanglement. In Figj3]we have presented the 
evolutions of junction B and D solving Eq. ( |39ll numerically. The left figure shows the 
evolution of junctions B and D both for matter and radiation dominated backgrounds. The 
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FIG. 4: Here the background evolution for F{x) and v{x) are presented with Fi = 0.5, = 
0.2,F'{x = 1) = 0.1. The solid (dashed) curves are for the radiation (matter) dominated back- 
grounds. 

right figure shows the length of the newly formed string fi^ stretching between junctions B 
and D which is S3 = S-f — . As mentioned previously for the junction to form we require 
that 5*3 > 0. Form the right figure we see that the junction is created in both cosmo logical 
backgrounds. After the junction formation, 5*3 reaches a maximum value indicating the 
unzipping of the newly formed strings BD. After this time, 5*3 reduces. When 5*3 = 
the junctions B and D meet again and the loops disentangle and pass by from each 
other. From the right figure we see that for the radiation dominated background the loops 
disentanglement indeed take place. However, for the matter dominated background, we see 
that before loops find the opportunity to disentangle, they shrink to zero. As explained 
below Eq. (1511) the junction B unzips sooner than junction D which is also demonstrated 
in the left figure of Fig. [31 

B. Loops with intermediate sizes 

For the loops with the sizes comparable to the Hubble radius at the time of collision, 
Fi ~ 1, we can only do numerical analysis. On the physical grounds one expects that the 
evolution of F{x) and v{x) is less sensitive to the background cosmological expansion as 
compared to large super-horizon size loops. In Fig. [4] we have presented the background 
evolution of F{x) and v{x) for Fi = 0.5. As expected v{x) changes slowly and F{x) evolves 
similarly for both matter and radiation dominated backgrounds. In Fig. [5]we have presented 
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FIG. 5: Here we present the evolutions of junctions B and D and the string BD for Fi = 0.5, a = 
7r/9 and k = 1.2. In the left figure the upper solid red curve represents whereas the lower solid 
blue curve is that of for the radiation dominated background. The dashed curves represent 
the corresponding curves in the matter dominated background. The right graph represents the 
length of the newly formed string BD, S3. The solid (dashed) curve is for the radiation (matter) 
dominated background. We see that in both backgrounds, after junction formation, the string 
BD reaches a maximum length and get unzipped. We also observe that the onsets of string BD 
unzipping and the loops disentanglement happen sooner in a radiation dominated era. 

the evolutions of junctions B, D and the string BD . For both radiation dominated and 
matter dominated backgrounds we see that the junctions are formed followed by the string 
BD unzipping and the loops disentanglement. One observes that the onsets of string BD 
unzipping (when 5*3 reaches a maximum) and also the loops disentanglement happen earlier 
for the radiation dominated era as compared to the matter dominated era. This may be 
interpreted by noting that for the radiation dominated backgrounds the loops reenter the 
horizon and shrink sooner as can be seen qualitatively from Eq. ( l42l) . 

C. Small sub-horizon size loops 

For small loops, Fi ^ 1, as explained before one expects the background cosmological 
evolutions do not play important roles. We have presented the analytical results for small 
loops in subsection IIII B[ In Fig. [6] we have presented the full numerical solutions of F{x) 
and v{x). As expected F{x) shows simple periodic behavior and v{x) does not change much 
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FIG. 6: Here the solutions of F{x) and v{x) for both matter dominated and radiation dominated 
backgrounds are shown for Fi = 0.01, fi = 0.2 and F'{x = 1) = 0. As expected, the background 
cosmological evolutions do not play important roles so F{x) indicates simple periodic behavior and 
v{x) changes slowly. 



in each period. In Fig. [7] we have presented the junctions evolution. As expected, the 
junctions evolutions are identical for bath matter and radiation dominated backgrounds. 
We also observe that the string BD unzipping and the loops disentanglement take place in 
this example. As demonstrated analytically in Eq. ( ISTl) the junction B unzips sooner than 
junction D which is also demonstrated in the left figure of Fig. [71 



V. CONCLUSION 



In this work we have studied the cosmic strings collision in cosmological backgrounds. 
After presenting the general formalism in section |Tl] we have concentrated to the example 



of colliding loops 
of simulation in 



he motivation for this work was to understand analytically the findings 
20| where it was found that there were little contributions from the bound 
states strings in their multiple strings network. One can understand this phenomena as fol- 
lows. For the junctions to develop upon strings collision, some appropriate initial conditions 
should be satisfied. These depends on the relative tensions of the colliding strings, the angle 
of collision and their relative velocities. Yet the more interesting observation is that even 
when junctions are created, they can not grow indefinitely and the bound state strings start 
to unzip. 
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FIG. 7: Here the junctions evolutions are shown for Fi = 0.01, q = 7r/9 and k = 1.2. In the left 
figure, the upper red (lower blue) curve shows S^{S^). Interestingly the curves corresponding to 
the radiation and the matter dominated backgrounds coincide to each other. As demonstrated by 
Eq. (|51|) the junction B unzip sooner than the junction D. The right figure shows the length of 
string BD given by 53. Again the curves corresponding to the matter dominated and the radiation 
dominated backgrounds coincide. 



As described in [19|, for straight cosmic strings at collision the junctions do not unzip 
once they are materialized. However, for colliding loops in a fiat background the zipping and 
unzipping generically happen The natural question is how sensitive are these results 

to the expansion of the Universe. Here we find some interesting results indicating that the 
background expansion plays important roles in strings zipping, unzipping and their eventual 
disentanglement . The key parameter here is the relative size of the loops compared to the 
Hubble expansion rate at the time of collision. For large super-horizon size loops one may 
approximate them with straight strings. This implies that if junctions are formed upon 
loops collision it will grow initially as in straight strings examples. However, as the Universe 
expands the loops stretch conformally until they re-enter the Horizon. Meanwhile their 
velocities reduces rapidly. The net effect is that the rate of bound state strings creation 
slows down until it starts to unzip. Eventually the loops disentangle from each other and 
pass by from each other in opposite directions if they did not shrink to zero by then. On the 
other hand, for small sub-horizon size loops one can neglect the effects of the expansion and 
the results of [3| holds true. The case of colliding loops with the sizes comparable to the 
Hubble radius at the time of collision is more non-trivial which shares some features with 
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small and large loops cases. 

Our numerical investigations also show that the junction formation and the zipping and 
unzipping phenomena are sensitive to the angle of collision a and the strings relative tensions 
parametrized by k. It would be interesting to study these phenomena in the multi parameter 
space of and a. Also to simplify the analysis, we have restricted ourselves to the 

example of coplanar colliding loops with equal tensions and radii. It would be interesting 
to study the general case where the loops have different sizes and orientations and there 
are hierarchies between the sizes of the loops and the Hubble expansion rate at the time of 
collision. 
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